Scattering Length Scaling Laws for Ultracold Three-Body Colhsions 
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We present a simple picture that provides the energy and scattering length dependence for all 
inelastic three-body collision rates in the ultracold regime for three-body systems with short range 
two-body interactions. In particular, we present the scaling laws for vibrational relaxation, three- 
body recombination, and collision-induced dissociation for systems that support s-wave two-body 
collisions. These systems include three identical bosons (BBB), two identical bosons (EBB'), and 
two identical fermions (FFF'). Our approach reproduces all previous results, predicts several others, 
and gives the general form of the scaling laws in all cases. 

PACS numbers: 34.10.-|-x,32.80.Cy,05.30. Jp 



The use of external magnetic fields to control the 
atomic interactions in trapped ultracold quantum gases 
has proven to be an extraordinary tool to explore dif- 
ferent quantum regimes. At low temperatures, only the 
two-body s-wave scattering length a is needed to charac- 
terize the atomic interactions and can assume practically 
any value from — oo to -l-oo by tuning a magnetic field 
near a diatomic Feshbach resonance. This tunability has 
been used, for instance, to convert fermionic atoms into 
weakly bound bosonic molecules which, in some cases, 
were remarkably long lived . In the quantum degener- 
ate regime, this system provides a unique opportunity to 
explore the crossover between a Bose-Einstein condensate 
(BEC) of molecules and the BEC of atomic Cooper pairs, 
the BEC-BCS (Bardeen-Cooper-Schrieffer) transition. 

Since the molecular lifetimes are infiuenced by three- 
body processes, these experiments underscore the im- 
portance of knowing the dependence of ultracold three- 
body collision rates on a. In particular, vibrational re- 
laxation, X+(XX)*^X+XX, releases enough kinetic en- 
ergy to free the collision products from typical traps, 
leading to molecular loss. Collision induced dissocia- 
tion, X-f XX-^X-fX-|-X, can also contribute to molecular 
loss for large a. The other possible three-body process, 
three-body recombination X-|-X-|-X^X-|-XX, contributes 
to atomic losses. In an experiment, however, the domi- 
nant process is determined in part by the threshold and 
scattering length scaling laws. 

Threshold laws, which give the energy dependence for 
small collision energies, dictate the partial wave that 
dominates for each process. The threshold law for elastic 
two-body collisions, for instance, gives an s-wave cross 
section that is constant at threshold while the ^-th par- 
tial wave is suppressed by a factor of E'^^ . Combined with 
permutation symmetry requirements, the threshold law 
thus leads to the familiar conclusion that the collision 
cross section for two indistinguishable fermions vanishes 
at low energies. Three-body threshold laws similarly de- 
pend on the number and kind of identical particles for 
each partial wave J'^ , where J is the total orbital angu- 
lar momentum and tt is the parity 

While general results exist for three-body threshold 



laws, no similarly general scattering length scaling laws 
have yet been obtained. Specific cases are known: the 
recombination rate for B-l-B-l-B collisions, where B is a 
bosonic atom, scales roughly as for 0+ and all a 0, 0| , 
and as for 2+ and a>0 5]. It has also been shown for 
a>0 that the relaxation rate for B-I-BB collisions is linear 
in a 0, but is a~^-^^ for F-I-FF' collisions 0, where F 
and F' are distinguishable fermionic atoms and a is the 
F-l-F' s-wave scattering length. This scaling law is qual- 
itatively different from the other known scaling laws and 
begs the question: how do the three-body processes for 
other systems scale with a? 

In this Letter we present a simple physical picture 
within which both the energy and scattering length de- 
pendence of all ultracold three-body collision rates (for 
short-range two-body interactions) can be derived and 
understood. A simple picture emerges because all three- 
body systems can be represented by one of four proto- 
type systems (two each for a>0 and a<0). In all cases, 
though, the rate limiting step is tunneling through a po- 
tential barrier in the initial channel. The barrier is deter- 
mined here from effective three-body potentials obtained 
from the adiabatic hyperspherical representation 0. A 
simple WKB approximation to the tunneling probability 
is then sufficient to give both the energy and scatter- 
ing length dependence. Besides their practical use, the 
scaling laws make explicit the pervasive infiuence of Efi- 
mov physics on ultracold three-body collisions. We 
restrict our discussion to systems with equal masses and 
to symmetries that support s-wave two-body collisions, 
representing most cases of experimental interest. 

In the adiabatic hyperspherical representation, the 
three-body effective potentials and couplings are deter- 
mined from the adiabatic equation 8], 



(1) 



where denotes all hyperangles and R is the hyperra- 
dius that, roughly speaking, gives the overall size of the 
system. The adiabatic Hamiltonian _ffad includes the hy- 
perangular kinetic energy as well as all interactions. By 
expanding the total wave function on the adiabatic basis 
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$1,, the Schrodinger equation is reduced to (atomic units 
will be used unless otherwise noted): 
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2^1 dR^ 
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In this expression, /i is the three-body reduced mass, E 
is the total energy, Fy is the hyperradial wave function, 
V^iji is the nonadiabatic coupling responsible for inelastic 
transitions, and W^, is the effective potential. 

For short range interactions, the asymptotic behavior 
of Wv can be derived analytically this limit, the 

molecular channels, which represent atom-molecule scat- 
tering, and the three-body continuum channels, which 
represent collisions of three free particles, are given by 



and Wi, — 
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respectively. The molecular bound state energy E^ii is 
labeled by the ro- vibrational quantum numbers v and I'; 
I is the relative angular momentum between the molecule 
and the atom; and A is a positive integer that labels the 
eigenstates of the hyperangular kinetic energy. 

When the s-wave scattering length \a\ is large, how- 
ever, an intriguing phenomenon known as the Efimov ef- 
fect occurs j9|, modifying the behavior of W^. These 
modifications must be considered in order to properly 
predict the dependence of the three-body rates on a 
0, 0, H, 0. In the hmit |a|3>ro, where rg is the char- 
acteristic size of the two-body potential, Eq. Q applies 
only for R^\a\. For ro^i?<C|a|, the potentials are pro- 
portional to -R^^, but can now be attractive as well as 
repulsive. Strictly speaking, the term "Efimov effect" ap- 
plies only to the emergence of an infinity of three-body 
bound states for |a|— >oo. We will instead use the term 
"Efimov physics" to indicate the qualitative change in 
behavior exhibited by any system whenever at least two 
of the three possible s-wave scattering lengths is large. 

Based on the modifications due to Efimov physics, 
we can classify all three-body systems into one of two 
categories: those with an attractive dipole potential for 
ro<C-R<tC|a| and those without. For equal mass systems, 
only 0+ bosonic systems fall into the first category. In 
these systems, the attractive potential appears in the 
highest vibrationally excited s-wave (or weakly bound) 
molecular channel (/'=0) for a>0, and in the lowest con- 
tinuum channel for a<0. The effective potentials for all 
higher channels are repulsive, but the coefhcients differ 
from those shown in Eq.® due to Efimov physics. These 
effective potentials are conveniently parametrized by, 
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The constants sq and s^ depend on the number of res- 
onant pairs as well as the number of identical particles. 
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So (Po) s„ (p„) 


BBB 


0+0 
1" 1 3 
2+2 2 


1.0062378 4.4652946 
2.8637994 6.4622044 
2.8233419 5.5082494 


EBB' 


0+0 
1" 1 1 


0.4136973 3.4509891 
2.2787413 3.6413035 


FFF' 


0+0 2 
1" 1 1 
2+2 2 


2.1662220 5.1273521 
1.7727243 4.3582493 
3.1049769 4.7954054 



TABLE I; Coefficients of the potentials in Eqs. Ex- 
cept for 0+ bosons, the constants correspond to po and p^. 



For all other cases, including J>0, the effective potentials 
are repulsive in the weakly bound and continuum chan- 
nels, for a>0 and a<0. Like the 0+ bosonic case, these 
potentials are parametrized by constants po and p,^: 
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In all cases, deeply bound molecular channels are essen- 
tially independent of a. 

The constants sq, s^/, po, and can be obtained an- 
alytically 1^ , and numerical values are shown in Table U 
for the two partial waves that dominate relaxation and 
recombination near threshold. Table also shows the 
minimum I and A [see Eq. allowed by permutation 
symmetry that, in turn, determine the dominant partial 
wave for the three-body rates near threshold 0- 

The relevant effective potentials and couplings are 
sketched in Figs. Q^a) and db) for 0+ bosonic systems 
and in Figs. Wic) and H^d) for all other cases. For 
ro<R<\a\, the potential curves are given by Eqs. Q and 
©; for i?>|a|, the curves are those for finite |a| in Eq. JSJ. 
The lowest continuum channel, labeled "a" , is the initial 
channel for recombination. Channel "/?" is the weakly 
bound molecular channel and is the initial channel for re- 
laxation. Channel "7" is a deeply bound molecular chan- 
nel in all cases. The nonadiabatic couplings are sketched 
in Fig. n and indicate the regions where inelastic tran- 
sitions are most likely. From our numerical calculations 
and on physical grounds, we believe that Fig.^represent 
all three-body systems near threshold. 

Knowing the dependence of the potentials and cou- 
plings on R and a enables us to derive the scaling laws 
from the definitions of the rates. 
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in terms of the transition probabilities ll/^p. Only the 
wave vector important for the threshold law has been in- 
cluded — P=2ii{E - Eyv) for Koi and k'^^2^iE for K3. 
The transitions proceed via tunneling in the initial po- 
tential to the R where the coupling peaks. The transition 
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FIG. 1: Schematic plot of the three-body effective potentials 
for 0^ bosonic systems for (a) a>0 and (b) o<0, and for the 
others systems, for (c) a>0 and (d) a<0. 



probability can thus be approximated by the WKB tun- 
neling probability (including the Langer correction 



iexp 



(7) 



The Langer correction is crucial for obtaining the correct 
scaling with E and a, and introduces "tunneling" even 
when there is no barrier in W^. At ultracold tempera- 
tures, the classical turning point is much larger than \a\. 

In all cases, the couplings peak at R^\a\, R^tq, or 
both. Relaxation for a>0, for instance, only occurs at 
small distances, i?«ro, where the coupling peaks [Vpj in 
Figs.QJa) and^c)]. It follows that the tunneling proba- 
bility in the initial channel (3 must be evaluated by inte- 
grating Eq. Q from the classical turning point rc to rg. 
This range, however, spans both kinds of potentials so 
that the tunneling in the two regions must be included. 
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Similarly, the transition probability for recombination 
must also include tunneling in two different regions be- 
tween Tc and tq. Recombination for a>0 can occur at 
two distances, though — R^ a and i? w rg — since the 
coupling peaks both places, leading to 
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where each term corresponds to a different reaction path- 
way. Although not indicated here, these paths can inter- 
fere and will be discussed below. For a<0, all couplings 
peak at small distances i?«ro. Applying the arguments 
above, relaxation and recombination are, respectively. 



|T,,P^p(fi,„ and |T,„P:.P(:1|„|P£1^,„. 



(10) 



We can now determine both the threshold laws and the 
scattering length scaling. Since Wi, has the same form 
in the region i?>|a| in all cases [Eq. |j2Jl], the tunneling 
probabilities in this region are also the same, 

« {kar^' and P^%^ cx (fca)^'+^ (H) 

Equation completely determines the threshold laws 
for each process The scaling with a, however, will 
be strongly modified by the tunneling probability in the 
region ro<R<\a\ due to Efimov physics. 

For 0+ bosons with a>0, the /3-channel potential is at- 
tractive in the region ro<R<a. Consequently, the prob- 
ability Pa^ro [Eqs. © and 0] is not a tunneling proba- 
bility, but rather a transmission probability that can be 
determined from general arguments based on the known 
solutions for Eq. I0J . This analysis gives 
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sinh(2?7) 
sin^ [so ln(a/ro)+<E>]-|-sinh^(r)) 



a, 



(12) 



since Z=0 gives the leading contribution in Eq. Hll|) . The 
constant rj is related to the probability for transitions at 
small distances P| and $ is an unknown small- i? phase. 
Equation (|12|l was deduced in Ref. 6], generalizing the 
result presented in Ref. by letting the transition prob- 
ability at small distances be less than unity. In both 
cases, however, the linear dependence on a was obtained 
by dimensional and physical arguments, while it follows 
directly from Eq. Hll|l in the present analysis. 

For 0+ bosons with a<0, the initial state for relaxation 
is a deeply bound vibrational state independent of a. The 
transmission probability — and thus the relaxation rate 
— does not depend on a either, leading to 



V^-ci OC A„k r. 



21 21+1 







(13) 



The constant will generically represent small- i? phy- 
sics that can give resonances due to a three-body Fesh- 
bach resonance, but is otherwise independent of a. 

For all cases other than 0"*" bosons, the initial channel 
for relaxation is never attractive in the region ro<R<\a\ 
[Figs.ffl^c) and^d)] so that relaxation proceeds by tun- 
neling only. For a>0, the relaxation rate is 



(14) 



The rate thus decreases with a whenever 2l+l<2po, 
yielding longer molecular lifetimes for larger a. For the 
specific case of a weakly bound FF' s-wave molecule, 
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J" 


E a > a < 


E a > a < 


BBB 


0+ 
1" 
2+ 


const a const 

^2 Q,"^'^^** const 
fe'' a-"''*^ const 


const (fc*) a* |a|* 


EBB' 


0+ 

1" 


const a const 

^2 a'^-'"^^ const 


const (fc*) a* |a|* 


FFF' 


0+ 

1- 

2+ 


const a~^'^^'^ const 

e a^^-s* const 
fc^ const 





TABLE II: Threshold and scattering length scaling laws for 
three-body rates. Boldface indicates dominant contributions. 



po=2.1662220 so that Vrci oc a"^-^^^'"^'*. This result agrees 
with the recent prediction of Petrov et al. and is con- 
sistent with experiments [IJ. Like the 0+ boson systems, 
relaxation does not depend on a for a<0, 
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The present analysis applies equally well to recombina- 
tion. For instance, recombination of O"*" bosons with a>0 
is determined from Eq. Q , which includes three different 
recombination paths. Only the interference between the 
first two terms in Eq. @ will be included here, since we 
expect the third term will be suppressed. In fact, this 
interference is well known 0, 0| , and the present analysis 
reproduces the known expression plus a modification: 
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The first term in Eq. (|16|l is the usual result, while the 
second is due to recombination at small R. Because the 
small- i? coupling lies so far into the classically forbidden 
region, we expect A^j to be small. For a<0, the present 
analysis yields the same expression found in Ref . jj] , 
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sinh(2?7) 
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(17) 



The recombination rate for all other cases can be deter- 
mined, giving for a>0 and a<0, respectively, 
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i+AA^-^r"+B,{^) 

2po 



2p„ 



\2X+A 



(18) 
(19) 
The 



predicting an asymmetry in i^a for a<0 and a>0. 
constants A^i and i?^ are expected to be small. 

Our numerical results for B+BB and F+FF' collisions 
confirm Eqs. ((T^ ~ ((T5|l . In fact, in our model for F-I-FF' 
the contributions for both s- and p-wave final states scale 
as a"^'^^, emphasizing that the scaling only depends on 
the initial state. Equations (|16|l - p8() have been verified 
numerically and their limitations studied Q . 



Tabled summarizes the scaling laws, showing only the 
main power-law behavior of each rate. The two dominant 
partial waves, determined by their energy dependence, 
are shown for each process. For completeness, we in- 
clude the scaling laws for the dissociation rate (for 
which k'^—2fiE). It is interesting to note that for F-I-FF' 
collisions where the relaxation rate decreases with a, the 
dissociation rate grows (see Table 1^) and will eventu- 
ally become the dominant mechanism of molecular loss. 
The table also indicates the competition from the next 
leading term for finite temperatures. For fermion relax- 
ation, this term is comparatively more important than 
for boson relaxation. 

In this Letter, we have deduced the scaling laws for ul- 
tracold three-body collision rates for all symmetries. We 
have used a simple and intuitive approach that describes 
all three-body collision processes in the same framework. 
For relaxation and recombination, the scaling laws de- 
pend only on the initial state. Any dependence on the 
final state is expected to enter via the coupling terms, 
and should at most be weakly dependent on energy and 
scattering length. The present results apply only in the 
threshold regime, i.e. when the collision energy is the 
smallest energy in the system j^. We have also shown 
the remarkable influence of Efimov physics — even for 
systems without bosons. 

This work was supported by the National Science 
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